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Abstract 

We study the associated production of charged Higgs bosons with W gauge bosons 
in high energy e + e~ collisions at the one loop level. We present the analytical results 
and give a detailed discussion for the total cross section predicted in the context of 
a general Two Higgs Doublet Model (THDM). 
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1. Introduction 



The discovery of the standard Higgs boson is one of the major goals of the present and 
future searches in particle physics. Direct or indirect results give stringent lower and upper 
bounds on its mass [I], [|. Moreover, the problematic scalar sector of the Standard Model 
(SM) can be enlarged and some simple extensions such as the minimal Two Higgs Doublet 
Model (THDM) versions || are intensively studied. The two most popular versions of 
the THDM (type I and II), differ in some Higgs couplings to fermions, but in both types 
and after electroweak symmetry breaking |7], the Higgs spectrum is the same. From the 
8 degrees of freedom initially present in the 2 Higgs doublets, 3 correspond to masses of 
the longitudinal gauge bosons and we are left with 5 degrees of freedom which manifest 
themselves as 5 physical Higgs particles (2 charged Higgs H ± , 2 CP-even H, h and one 
CP-odd A). The charged Higgs H ± , because of its electrical charge, is noticeably different 
from the other SM or THDM Higgs particles. In other words the discovery of a charged 
scalar Higgs boson would attest definitively that the Standard Model is overcome. On the 
contrary a neutral Higgs experimental evidence has to be refined to go beyond the SM[Q ||. 

Furthermore, one can construct models with an even larger scalar sector than in the 
THDM, one of the most popular being the Higgs Triplet Model (HTM) [Q. A noteworthy 
difference between THDM and HTM is that the HTM contains a tree level ZW^H^ 
coupling while in the THDM this coupling is only generated at one loop level [[3], [|. In 
models of the HTM type, the physical predictions are sometimes very different from typical 
ones in the SM and or the THDM. In our study, we will only consider models with two 
Higgs doublets. 

The study of THDM is also very interesting insofar as it provides a general framework to 
test the minimal suspersymmetric extension of the SM (MSSM) ||. Of course the particle 
spectrum of the MSSM doubles that of the THDM, but the scalar Higgs sectors, responsible 
for the electroweak symmetry breaking, are analogous, up to some mass constraints and 
couplings in the MSSM. In any case, the study of the various production mechanisms of 
charged Higgs bosons would give information about new physics beyond the SM. 

In future colliders or linear accelerators devoted to the electron - positron annihilations, 
the simplest way to get a charged Higgs is through H ± pair production. Such studies have 
been already undertaken at tree-level and one-loop orders |TIJ [T| and shown that e + e~ 
machines will offer a clean environment and in that sense a higher mass reach, especially 
if the 1-2 TeV options are available. We mention also that charged Higgs bosons pair 
production through laser back-scattered 77 collisions has been studied in the literature 
fl"2|| and found to be prominent to discover the charged Higgs boson. 

Note also that charged Higgs bosons can be produced in hadronic machines: 

(i) single charged Higgs production via gb — > tH~ , gg — > tbH~, qb — > q'bH~ |13 



(ii) single charged Higgs production in association with W gauge boson via gg — > 
W ± H^ or bb -> W ± H^ @ 

(Hi) H^ 1 pair production through qq annihilation would be feasible only if the charged 
Higgs decay in top-bottom is kinematically forbidden |15 |, otherwise one would have to 
look at bosonic decay modes (ex. — > h°W , H ± — > A !-^, etc.). 
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The aim of this paper is the study of the production of the charged Higgs boson in 
association with a W gauge boson through electron positron annihilation. This process is 
kinematically better suited than the if ± pair production when the Higgs mass is larger 
than the W mass. Although it is a rare process in the THDM, loop or/and threshold 
effects can give a substantial enhancement. Moreover, once worked out, any experimental 
deviation from the results within such a model should bring some fruitful information on 
the new physics. 

The paper is organized as follows. In section II, we review all the charged Higgs boson 
interactions (with gauge bosons, scalar bosons and fermions), section III contains the 
notations and conventions while section IV is devoted to the on-shell renormalization 
scheme we will use. In section V we present our numerical results, and section VI contains 
our conclusions. 



2. Charged Higgs boson interactions 

2.1 Charged Higgs boson interactions with scalar bosons 



In this section we will follow the notation of ref |16 . Recalling the most general (dimension 



4) SU{2) wea k x U{1)y gauge-invariant and CP-invariant scalar potential for the THDM 
V($i, $ 2 ) = Ai($+$i - v 2 x f + A 2 ($+$ 2 - v\f + A 3 ((^$i - v\) + ($+$ 2 - v\)) 2 
+A 4 (($^i)(^$2) - ($^$ 2 )(<I> 2 f $i)) + A 5 (i?e($+<l> 2 ) - Vl v 2 ) 2 + 
A 6 [/m($+$ 2 )] 2 + A 7 , (2.1) 

we have the doublet fields $i and $ 2 with weak hypercharge Y=l, the corresponding 
vacuum expectation values v\ and t> 2 , and the coefficients Aj as real- valued parameters. We 
will assume the arbitrary additive constant A7 to be vanishing. Via the Higgs mechanism, 
the W and Z gauge bosons acquire masses given by m 2 ^ = \g 2 v 2 and m| = \{g 2 + g' 2 )v 2 , 
where g and g' are the SU{2) wea k and U(l)y gauge couplings and v 2 = v \ + v\. The 
combination v\ + v\ is thus fixed by the electroweak scale through v\ + v\ = (2y / 2G'^)~ 1 , 
and we are left with 7 free parameters in eq. (|2.1|) , namely Aj and tan (3 = f 2 /fi. As shown in 



TT| , using straight-forward algebra, one can relate the coefficients Aj of the scalar potential 
to the masses of the physical Higgs bosons h, H, A, [|l6l in the following way: 



g 2 m 2 H ± g 2 sin2a 2 2 g 2 m\ 

A4 = ^4T ' As = 2<:im^ (m --^ ) " 4X3 ' Ae= 2<: (2 - 4) 
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Consequently one can choose as free parameters of the Higgs sector: the four Higgs boson 
masses (m H ±, uih, ?nh, ttia), «, tan/3 and A3. 

The trilinear self couplings required for our study are listed in the appendix A, Eqs. ( |A.1| )- 
( |A.7| ). The trilinear vertices H°H + H~ and h°H + H~ depend, besides on the Higgs boson 
masses, on A3 and tan/3. From Appendix A, eqs.( |A.l| , |A.2j ), it is obvious that one would 
get into conflict with the requirements from perturbative unitarity when tan/5 and/or A3 
are large. In our analysis we will take into account the following constraints when the 
masses and the coupling parameters are varied: 



Unitarity constraints will be respected in a simplified way, following |T7|], by imposing 
the condition 

\HHH\ < J^-(lTeV) 2 (2.5) 
1 1 2m w K J K ' 

on each trilinear scalar vertex HHH entering the process at one loop. 

Lower bounds on tan/3 in the general THDM have been obtained from the experi- 
mental limits on the processes e + e~ — > Z* — > h 0r y and/or e + e~ —>■ A 0/ ~f ||18|| . For light 
h masses, these bounds can be rather low jlSj . For our study we will restrict the 
discussion to values tan/5 > 0.5. 



The extra contribution 8p to the p parameter [H5| should not exceed the current 
limits from precision measurements |2T 



-0.0017 < 5p< 0.0027 



2.2 Charged Higgs boson interactions with fermions 

In the two-Higgs-doublet extension of the Standard Model there are different ways to couple 
the Higgs fields to the fermions. Conventionally they are classified in terms of the following 
categories, labeled as type I and type II models: 

i. Model type I: All quarks and leptons couple exclusively to the second Higgs doublet 
$2, with the coupling structure copied from the Standard Model. $2 gives mass to 
both up- and down-type quarks, invoking the charge-conjugate of $2 for the up-type 
quarks. Since supersymmetry forbids the appearance of the complex conjugate, type 
I models cannot be realized within the MSSM. 

ii. Model type II: To avoid the problem of flavor changing neutral currents (FCNC) PT| , 
one assumes that $1 couples only to down-type quarks and charged leptons and $2 
to up-type quarks (and eventually to neutrinos). The type II model is the pattern 
found in the MSSM. 

The general structure of the charged-Higgs boson interaction with a doublet of up- and 
down-type fermions is given by the vertex 

\H-„4 = J^_ {ja Uz2!) + + . (2.6) 
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In the specific models mentioned above, the Yukawa couplings have the form 



Yh = — ^ and Y3 = m d tanf3 for model II . (2.7) 
tan p 

V u d is the CKM matrix element, which we will approximate by unity. Models I and II lead 
to similar results as far as the effects originating from the top-bottom loops for tan j3 close 
to 1 are considered, owing essentially to the form of the Yukawa couplings Y R and Y L . 



2.3 Charged Higgs-boson interactions with gauge bosons 

The interactions between the charged Higgs bosons of a two-doublet model and the gauge 
bosons are completely dictated by local gauge invariance and thus independent of the 
assumption whether the model is supersymmetric or not. These interactions follow from 
the kinetic term, involving the covariant derivative, in the Higgs Lagrangian 

£(i?„s,) + (i>„s<) = 

i 

+ igT a W- + ig'^LB^+Kd, + igf a W- + ig'^B,)^} , (2.8) 

where T a is the isospin operator and Y$ i the hypercharge of the Higgs fields; W a ^ denote 
the SU{2)l gauge fields, the U(l)y gauge field, and g (resp. g') the associated coupling 
constants. The individual couplings in the various vertices, derived from eq. ( |2.8| ), are 
listed in Appendix C. The , yH Zf W ± and the ZH^W^ vertices vanish at the tree level. 
Therefore, tree-level contributions to e + e~ — > W^H^ come only from neutrino exchange 
in the t-channel and from CP-even Higgs mediated s-channel. All these contributions 
are strongly suppressed by the small electron mass. Hence, the process under study is 
essentially a loop-mediated process. 



3. Notations and conventions 

We will use the following notations and conventions. The momenta of the incoming electron 
and positron, outgoing gauge boson W + and outgoing Higgs boson H~ are denoted by pi, 
P2, k\ and k 2 , respectively. Neglecting the electron mass, the momenta, in the e + e~ center- 
of-mass system, are given by: 

Pl,2 
h,2 



2 



1,0,0,±1) 



1 ± ^ m K ±k sin 9, 0, ±k cos . 

2 v s 
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where \fs denotes the center of mass energy, 9 the scattering angle between e and W + , 
and k is determined by 

k 2 = (s — (m H ± + m w ) 2 )(s — (m H ± — m w ) 2 )/s 2 . 

The Mandelstam variables are defined as follows: 

s = (pi + P2) 2 = (h + k 2 f 

■ s s 

t = ( Pl - kx) 2 = (pa - k 2 ) 2 = -{jn 2 w + m 2 H ±) - - + -kcosO 

\ s s 

u=(pi- k 2 ) 2 = (P2 - h) 2 = ^(mw + m H±) ~ 2~ 2 KCOS ° 

s + t + u = mly + m 2 H ± . 



The differential cross section reads, expressed in terms of the one-loop amplitude M 1 : 



'>"'e+e-^W ± HT) = 7 JL n Y,\Mr- (3.1) 

Pol 



Ai 1 can be decomposed into six invariant matrix elements Ai and scalar coefficient func- 
tions Mi as follows: 

M^j^MiA. (3.2) 

i=i 

The basic matrix elements Ai are given by the following expressions: 

A = v(p 2 ) 4(h) 1 ^ 75 m(pi) 
A = ^(p 2 ) A*i)^^«(pi) 

As = v(p 2 ) fa 1 Y^ u(pi)(pi.e(fci)) 
1 — 7s 

A = ^(p 2 ) A - 2 u (Pl)(Pl- e ( fc l)) 

1+75 

A = u(P2) A 2 «(pi)(p 2 .e(A;i)) 

A = u(P2) A - 2 — w (Pi)(P2.e(fci)) . (3.3) 

The squared amplitude, after summation over the polarizations of the W + boson, gets the 
following form: 

^ - 2 -~tu) (A]KA |2 , Ai 



|2 



E l-^T = 2 S (|^i| 2 + |-M 2 | 2 ) - (m ^ ±m ^ tM) {4|A^ 1 | 2 + A\M, 

Pol Am W 

+4(m 2 w - t)Re[MiMl + M 2 M\] + (m 2 w - t) 2 [\M 3 \ 2 + |.M 4 | 2 ] 
+4(m 2 w - u)ite[A<iA<S + A^A^] + (™ 2 W - u) 2 [\M 5 \ 2 + \M 6 \ 2 ] 
-2(m 2 H± m 2 w + m 2 w s - tu)Re[M 3 M* 5 + M 4 M* 6 }} (3.4) 



The coefficient functions Aii can be read off from the explicit expressions in Appendix C, 
according to the decomposition ( |3.2|) . 



4. On-shell renormalization 

We have evaluated the one-loop induced process e + e~ — > W ± H T in the 't Hooft - Feynman 



gauge, and using dimensional regularization [22]. The types of Feynman diagrams are 
depicted in Figure 1. It displays the corrections to the ■j-W^H^ and to the Z-W^H^ 
vertices, box diagrams, contributions coming from the various mixings H~-W + and H~ 
G + in t and s channels, and finally the counter-terms. Note that the s-channel H~-W + 
mixing (Fig. 1.16 ) vanishes for an on-shell transverse W gauge boson. There is no 
contribution from the W + -G~ mixing because the j-G ± H T and Z-G ± H T vertices are 
absent at the tree level. All the Feynman diagrams have been generated and computed 
using FeynArts |^3] and FeynCalc |23| packages. We also use the fortran FF-package |25[ 



in the numerical analysis. 

Although the amplitude for our process vanishes at the tree level, complications like 
tadpole contributions and vector boson-scalar mixings require a careful treatment of renor- 



malization. We adopt the on-shell renormalization scheme of [26], for the Higgs sector, 



which is an extension of the on-shell scheme in f27j] . In this scheme, field renormalization 
is performed in the manifest-symmetric version of the Lagrangian. A field renormaliza- 
tion constant Z§ 12 is assigned to each Higgs doublet <3>i 2 . The Higgs fields and vacuum 
expectation values Vi are renormalized as follows: 

v^iZ^^-dv,). (4.1) 



With these substitutions in the Lagrangian ( |2.8| ), expanding the renormalization constants 
Zi = 1 + SZi to the one-loop order, we obtain all the counter-terms relevant for our process: 
a counter-term for the W — H 2-point function, and counter-terms for the 'yWH and ZWH 
vertices, visualized in Fig.1.26 to Fig. 1.28 (ki denotes the momentum of the outgoing W + ), 

6[W±H*\ = - hi A (4.2) 
5[A„W±H*\ = ieg^A (4.3) 

6[Z V W±H*) = ieg» v —A (4.4) 
c w 

with 

A = ——m w [ + -(8Z 02 -dZ* x )] . (4.5) 

In the on-shell scheme, the following renormalization conditions are imposed: 

• The renormalized tadpoles, i.e. the sum of tadpole diagrams T^h and tadpole counter- 
terms 5h,H vanish: 

T h + 5t h = 0, Tjj + Stjj — . 
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These conditions guarantee that vi j2 appearing in the renormalized Lagrangian Cn 
are located at the minimum of the one-loop potential. 

• The real part of the renormalized non-diagonal self-energy E#-w+(/c 2 ) vanishes for 
an on-shell charged Higgs boson: 

Re£ H -w+(m 2 H ±) = (4.6) 

(the * labels the renormalized quantity which is given in Appendix C). This renor- 
malization condition determines the term A, and consequently ^f^W^-fT* 1 ] and 
b\Z v W^H^\ are also fixed. The explicit expressions are given in Appendix C. 

The last renormalization condition is sufficient to discard the real part of the H~-G + 
mixing contribution as well. Indeed, using the Slavnov-Taylor identity 



k 2 Y> H -w+(k 2 ) — mwEH-G+ik 2 ) = at k 2 = m 2 H ± 
valid also for the renormalized quantities, together with eq. Q4.6|) , it follows that 

ReL H - G +(m 2 H ±) = . 

In particular, the Feynman diagrams depicted in Fig. 1.25 will not contribute with the 
above renormalization conditions, being purely real valued. 

The amplitudes corresponding to the two counter-terms (photon and Z boson ex- 
changes) in eqs. (fD| |4.4j) contain only the matrix elements A\ and A 2 as follows: 



Ml T = --A(A 1 + A 2 ) 

2 

Ml T = 6 — S -^A((g A + g v )A 2 + (g v - g A )A 1 ) (4.7) 

s — m z z cw 

with the electron-Z coupling constants 

^ V 4s w c w ^ w) i 9a ^ SwCw ^ 

and sw = sin 9w, cw = cos 9w- As a consequence, only the form factors M.\ and M. 2 can 
contain UV divergences; each of the remaining form factors should be UV finite by itself, 
which provides a good analytical check of the calculations. 



5. Numerical results and discussion 

The difference between the predicted cross sections in Model-I and Model-II is essentially 
due to the fermion loops containing the different Yukawa couplings. For small values of 
tan f3 the cross sections are approximately the same for both models of type I and II (Figure 
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2a). For tan/5 > 10 (Figure 2b), the cross section decreases to very small values in Model- 
I, corresponding to the Yukawa coupling proportional to 1/ tan/?. In Model-II, the cross 
section has a minimum for tan j3 ~ 30 and increases again for large tan /3, following the b 
Yukawa coupling proportional to tan/?. 

In the following more explicit discussion we consider the range 0.5 < tan/3 < 10, where 
both type I and II models give practically the same numerical results. In this interval 
for tan/? the total cross section varies from 7.5 fb to 0.012 fb, for ^/s = 500GeV and 
m H ± = 220GeV. 

We will take the following experimental input for the physical parameters [3(| : 



2 

• the fine structure constant: a = |- = 1/137.03598. 

• the gauge boson masses: Mz = 91.187 GeV and Mw = 80.41 GeV. In the on-shell 
scheme we define sin 2 ^ by sin 2 6^ = 1 — so that it receives by definition no 
loop corrections. 

• the top and bottom quarks masses are taken to be: m t = 175 GeV and = 4.5 
GeV. 

Fig. 3 shows the top-bottom contribution to the integrated cross section as a function of 
the center of mass energy yfs for four values of mn± and for two values of tan/3, 0.6 and 
1.6. It can be seen that for a small tan/3 (Fig. 3. a) the top quark effect is enhanced. One 
can reach a cross section of 3.5 fb for a charged Higgs mass of about 220 GeV. The cross 
section is enhanced for m#± = 170 GeV owing to the proximity of the normal threshold 
cut of the three-point function at rrtn± = m t + rrib. 

As tan j3 increases the top quark effect decreases, leading to almost an order of mag- 
nitude suppression of the cross section for tan/3 = 1.6 (Fig. 3.b). For large values of 
tan /3(~ 50) the bottom quark contribution becomes leading and of comparable magnitude 
to that of the top quark in the small tan /3 region. 

For the general THDM we will present our numerical results in the following three 
configurations (where all the masses are in GeV): 

• C\\ m H ± = 220, rriH = 180, = 90, rriA = 80, tana = 1.4, tan/3 = 3.6 

• C 2 : m H ± = 300, m H = 280, m h = 120, m A = 220, tana = 2.4, tan/3 = 1.6 

• C3: m H ± = 400, ran = 380, = 120, uia = 370, tana = 3, tan/3 = 2 



In those three cases the present experimental bound on 5p is respected. In Fig.4.a,b,c 
we show the total cross section (including all virtual boson and fermion contributions) as 
a function of the center of mass energy yfs for the three configurations listed above and 
for four different values of A3 which satisfy the nominal unitarity constraints ( |2.5| ) One can 
see from those curves that, for a fixed tan/3, the cross section increases with A3. The effect 
of A3 arises essentially via the trilinear vertices H°H + H~ and h°H + H~ . Note that values 
of A3 larger than the ones chosen in Fig.4.a,b,c start violating the unitarity constraints on 
the H°H + H~ and h°H + H~ couplings. 
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Fig.4.d shows the (A3 independent) box contributions to the cross section in the Ci j2 ,3 
configurations. This contribution remains generically suppressed (0.01-0.005 fb for s/s =500- 
1000 GeV), even in the favourable light Higgs mass configurations such as in case C\. 



Figures 5. a, 5.b and 5.c show the total cross section as a function of tan/3 in the Ci.2,3 
configurations with A3 = 0.1, and for yfs = 500 Gev, 1 TeV and 1.5 TeV respectively. One 
can see that the cross section reaches a minimum for moderate values of tan f3 while it gets 
much larger (a few orders of magnitude) for both low and high values of this parameter. 
The first enhancement (small tan ff) is due to the top quark effect as discussed above, while 
the second enhancement comes from the effect of large tan (3 in H°H + H~ and h°H + H~ 
couplings. With the parameter choice of Fig.5, the nominal unitarity condition (|2.5|) forbids 
tan /3 larger than 14. One should, however, stress that a different choice of the Higgs masses 
and A 3 closer to the MSSM configurations would suppress the Higgs sector effect in the 
large tan/3 regime, thus allowing even larger tan/? values to be consistent with Eq. (|2.5| ). 
In this case large effects from the bottom quark sector become dominant. 

In Fig. 6 we show the dependence of the cross section on the charged Higgs mass for 
y/s = 500 Gev, 1 TeV and 1.5 TeV in the case where tan/3 = 2, tana = 3, A3 = 0.1, 
mil = 180 GeV, rrih = 90 and uia = 80 GeV. One notes the high sensitivity of the 
cross section to the Higgs mass when the latter is above the W~h and hi decay thresholds 
(which happen to roughly coincide owing to our input values, and correspond to the kink 
in the plots). Below this threshold, the cross-section is much less sensitive to the Higgs 
mass, however, one should keep in mind that this region corresponds to the opening of the 
charged Higgs pair production which yields a much more interesting event rate. 

We mention that our results are in qualitative agreement with those of ref.|32|], inasmuch 
as we can compare. Note however that we did not assume any Higgs mass rules, and also 



the renormalization scheme in [32] was not explicitly defined. Quite recently, after this 



work was completed, another paper |33|] appeared dealing with the same subject. 



6. Conclusions 

To summarize, we have computed the associated production of a charged Higgs boson 
with a gauge boson W at high energy e + e~ collisions. The calculation is performed within 
dimensional regularization in the on shell scheme. The study is done in the general two 
Higgs doublet model taking into account constraints on the p parameter and unitarity 
constraints on the trilinear Higgs couplings. We have shown that in the small tan (3 regime 
the top effect is enhanced leading to important cross sections (about 1 fb), while the 
leading contributions in the large tan/3 regime, come from the trilinear HH ± H T and 
hH ± H :f vertices. If the charged Higgs is too heavy to be produced in pairs in e + e~ future 
machines, the associated production with a W boson would be the only means to look for 
direct evidence for it. The smallness of the cross section would require, however, a very 
high luminosity option. 
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Appendix A: 

In this section we list the Feynman rules for the 3-point vertices involving charged Higgs 
bosons, in the general THDM. 

• 9 r in w 2 m \\ S i n ( a + (3)^A^ 2 1/2 2\ 

9hoh + h- = -t-Z-[cos(j3 - a )(m 2 H± -^f ) + ^—f±{4\ 3 v 2 + - (m 2 H + m 2 h ) 
m w 2 sin2p 2 

' 'sin2a + 2sin(a - 0) cos(a + (3))(m 2 H - ml)}) (A.l) 



2 sin 2/3 

• # r • ,n w 2 COS ( a + ^)r,A 2 1/2 2\ 

TTljy Z Sill Lp Z 

' 'sin 2a + 2sin(a + /3) cos(a - l3)){m 2 H - m 2 h )}} (A.2) 



2 sin 2/5 

-i£ sin(/3 - a) (m 2 H± - m 2 H ) 
9h°h±g* = o l A -3j 

igcos(/3-a)(m? H± - m|) 

9h°H±G* = o (A 4 ) 

2m w 

m?r ± — m 2 A , . _ N 

0aoh±gt = =F — — /= — (A.5) 



2 cos(/3 — a) 



9h°g+g~ = -im H -= (A.6) 

V2w 



2 sin(/3 — a) 



9hPG+G- = -im h — -j=^ — (A. 7) 

The parameter A3 enters only gn°H+H- and gh°H+H-, while the vertices 9b°h±gt , gh°H±GT 
and (7yio#± G =F have a particularly simple form, proportional to m 2 H and m|, respectively. 

Appendix B: One-loop functions 



Let us briefly recall the definitions of scalar and tensor integrals |3T| we use. The inverse 
of the propagators are denoted by 

d = q 2 - ml , dt = (q + Vif ~ m 2 
where the Pi are the momenta of the external particles (always incoming). 
One-point function: 

Mm ^ = w J d q T 

where p is an arbitrary renormalization scale and D is the space-time dimension. 
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Two-point functions: 



Qk^V /WW 



f d D q±—^ 
J d§d\ 

Using Lorentz covariance, one gets for the vector integral 

with the scalar function Bi(p\,ml,m\). 
Three-point functions: 

J (IqUi-i 



(2tt/x)( 4 d > r d {hq^q^} 
°0, M ,^ - —2 / a 1 771 

111 J (ln(li(l<i 



'•2 



where pfj = (pi +pj) 2 . Lorentz covariance yields the decomposition 

C lt = Pl/t C 1 + P 2 li C 2 (B.l) 

C^ u = g^uCoo + PluPXvCn + P2fiP2„C 2 2 + {V\ixV2 v + P2nPl v )Cvi (B.2) 

with the scalar functions Ci^{p\ ,pf 2 ,p 2 , m o, m i, 
Four-point functions: 

V0,WV ~ — / « g A A A A ( B - 6 > 

Again, Lorentz covariance allows the decomposition 

Dfj, = pi M Di + p 2/ ,D 2 + p 3M L> 3 (B.4) 
= g^Doo + VXfxViv'Dii + P2 l xP2uD 22 + P^PzvDm + (pitf>2v + P2 l iP\u)D 12 + 

(Pl^PSu + P^Plv)D Vi + (p 3 ^P 2 u + P2»P3v)D 23 . (B.5) 

with the scalar 4-point functions Diij(j^,p\ 2 ,p\ z ,pl,p\,p\ 3 ,ml,m\,m\,m^). 

All the tensor coefficients can be algebraically reduced to the basic scalar integrals Aq, 
Bq, Cq and Dq. The analytical expression of all the scalar functions can be found in pi|, 
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Appendix C: One- loop amplitude 

In this appendix, we use the fermion-vector boson coupling constants as defined in terms 
of the neutral-current and charged- current vertices 

fL 1-75 , /fll + 75. 



[ffv»] = + g{r-7^) 

1 

\[2s w 



in the following notation: 



& = 9? = ~e, 

g uL = (1 - 2s 2 w e u ) = 2s 2 w e f 

4swCw 4swCw 



gf 



2s 2 w e d ) 



4s w c w 
2 

-2e d = -\e\, e e = -\e\. (C.2) 



Vertex diagrams 
Fermionic loops 
d-d-u exchange 

The diagram with the ddu triangle, Fig. 1.1, yields the following contribution to the one- 
loop amplitude Q3.2p for each V boson exchange (V — 7, Z): 

M1.1 = . J Nca ], z ( {~{m u Y u L d gf + m d Y R (g d v L - g d v R ))B (s, m 2 d , m 2 d ) - 

m u (m 2 u Y u L d g dL + m d m u Y* (g d v L - g dR ) - m 2 d Y u L d g dR )C + 
2(m u Y u L d + m d Y^)g dL C 00 }[gf A, + gf A 2 ) 

-{(m u (m 2 w + u) Y u L d g dL + m d Y* (m 2 w g dL - t g v R ))d + m u u Y u L d g d v L C 
+ (m u (m 2 H± + u) Y u L d g dL + m d Y* (u g dL - m 2 H± g dR ))C 2 } gfA x + 
{K (m 2 w + t) Y u L d g dL + m d Y* (m 2 w g dL - u gf )]d - m u t Y u L d g dL C + 
K (~<± ~ t) Y u L d gf + m d Y u R d (m 2 H± gf - t gf)]C 2 } gf A 2 
-2{gf (m u Y u L d C + m u Y u L d C x + (2m u Y u L d + m d Y R )C 2 + 
(m u Y u L d + m d Y R )C 12 + (m, Y u L d + m d Y R )C 22 )}[gfA 3 + gfA 6 ] 
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+2{m d Y* gf C x - g d v L {m u Y u L d C 2 + 

(m u Y u L d + m d Y u R d )(C 12 + C 22 ))}[#A + gfAs] ) . (C.3) 

Therein, all the Cf and Cy have the same arguments: (m 2 ^, s,m 2 H± ,m 2 ,m d ,d d ). Sum- 
mation has to be performed over all fermion families; in practice only the third quark 
generation is relevant. 



u-u-d exchange 



The corresponding expression for the the diagram with the uud triangle in Fig. 1.2 is ob- 
tained from the previous one by making the following replacements: 



Y 



L 



- ud 

and also A 



Y, 



R 

ud ' 



rrin <- 



An and A 



m d , 
-> A- 



uL 

9v 



dR 

9v 



9v 



l R and t 



-» u 



Bosonic Loops 



Listed are always the analytic expression for each generic diagram of Fig.l. The sum 
over all the particle contents, given in the associated tables, yields the corresponding con- 
tribution to the one- loop the matrix element (|3.2|). 



Fig.1.3 



Ml 



« 9vw+w 



m 



v 



2 9w+w-s9h-w+s{ ( -B (s, m\, m 2 w ] - (m 2 s - m 2 w + s)C + 



[m 



H± 



m 



w 



s)Ci + {m 2 H± 



m 



w 



+ s)C 2 + C 00 ) [ g fA 1 +g e v L A 2 



+(4d + C 2 - C l2 - C 22 )[gf {A 3 + As) + g e v L (A + A)]}- (C.4) 

All the Cj and have the same arguments: (m 2 ^, s, m% ± , m|, m 2 ^, m^). The couplings 
are summarized in the following table; S is a generic notation for one of the Higgs bosons. 



S 



h 



H 



9zzs 



m z 



s 0a 



m z 



9w+w~s 



m w 



m w 



Cf3a 



9h-w+s 



^0a 

2sw 



2s v 



9zw+w- 



s w 



9-yW+W- 



-1 



Therein, sp a = sin(/3 — a) and cp a = cos(/3 — a) 



Fig.1.4 

M 1A 



2a 2 



-9vSiSj9w-w+Sj9H-w+Si {Coo[9v Ai + 9v A 2 



s — mi 



-(2C + 2d + 3C 2 + C 12 + C 22 )[gf (A + As) + 9v (A + A)]} (C.5) 



eX 
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where the arguments of Cj and Cy are now as follows: (rri^, s,m^ ± ,m^,mg.,rrig.). The 
couplings are summarized in the following table: 



(Si,Sj ) 


gzSiSj 


9H-SiW 


gw-w+Sj 


(A Q , h) 


c f3a 


1 

2s w 




(A), H) 


Spot 
2sw c W 


1 

2s\y 





Fig.1.5 

a 2 

Afl.5 = 29VS i S j 9w+S j S k 9H-S i S k {Coo\gv t ^1 + # -A 2 ] ~ 

(C 2 + Ci 2 + C 22 )[# (A 3 + A 5 ) + # (A 6 + A)]} (C.6) 

All the Ci and have the same arguments: (m^, s, m%±, m 2 Sk , m|., m| ). The couplings 
are summarized in the following table: 



{Si,Sj,Sk ) 




gzSiSj 


9H±S t S k 


9w+SjS k 


(A , /i, G ± ) 







Wh-a g+ 


Spa 

2s\y 


(A , if, G ± ) 





i Sf3a 

2s w c w 


Wh+a g- 


Cf)a 

2sw 


(G ± , G ± , /i ) 


-1 


cos(20vk) 
2swCw 


9H-hG+ 


Spa 
2s\y 


(G ± , G ± , 5- ) 


-1 


cos(2^vi/) 
2syi/cw/ 


9H-HG+ 


c /3a 

2sw 


(ff±, ff±, ) 


-1 


cos(20via) 

2syi/Cyi/ 


9H-hH+ 


c /3a 
2s\y 


(ff±, ff±, ff ) 


-1 


cos(2^vk) 
2sv^cv^ 


9H-HH+ 


Sf)a 

2sw 



where the gu±SiS k couplings have been defined in appendix A. 
Fig.1.6 

a 2 

M L6 = 2 -9vv's i gw+v'S j gH-s i s j C [g e v R A^ + gf A 2 ] (C.7) 

s — m v 

where Co = Co(m^, s, m 2 H± , m 2 Si1 m 2 s .,myi). The couplings are summarized in the following 
table: 



{Sifty ) 


9ivs % 


9zv'Si 


9H+S.S, 


gw+vsj 


(K G ± , z ) 





m z ^^- 

* s w c w 


gn-G+h 


-m z s w 


(ff, G ± , Z ) 





c 0a 

m z 


gu-G+n 


-m z s w 


(G ± , h, W ± ) 


m w 


—m z s w 


gn-G+h 


m w 
S P« s w 


(G ± , ff , W ± ) 


mw 


-m z s w 


gu-G+n 


C P a ~s~w~ 
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Fig.1.7 



Mi.. 



2a 2 



-9vv's i gw+s % s ] 9H-v'sA~CoQ[g'^ i A\ + g e y A 2 \ + 



s — m 



v 



-C 2 + C 12 + C 22 )[gf (A 3 + As) + 9v (A + A)]} 



eL 



(C. 



All the Cj and Cy have the same arguments: (m^, s, ra# ± , ra|., ra|., ra V /). The couplings 



are summarized in the following table: 



(S^Sj ) 




9zw+Si 


gn-w+Sj 


gW+SiSj 


(G* fc) 


m w 


-m z s w 


2sw 


Spa 

2sw 


(G* H) 


mw 


—m z s w 


«/3a 
2sw 


c pa 

2sw 



Fig.1.8 



a 1 



rrir 



-9zzs9zw-w+9H-w+s{(B (s,m 2 s ,m 2 z ) + (2m 2 H± + 3m 2 w - 2s)C 



+ [m 2 H± + 3m 2 w -s]d + [3m 2 H± + m 2 w - s]C 2 - C m )[g e v R Ai + g e v L A 2 ] 
(2C - 2d + 3C 2 + C 12 + C 22 )[g e v R (A 3 + A 5 ) + gf (A G + A 4 )]} (C.9) 

where all the Cj and have the same arguments: (m^, s,m 2 j ± ,m 2 v ,m 2 z ,m'g). The 
couplings are summarized in the following table: 



s 


gzzs 


9h-w+s 


9zw+w- 


h 


SQa 

swcw 


<^pa 

2sw 


sw 


H 


* s w c w 


-Spa 

2sw 


sw 



Fig.1.9 

Mi.g 

Fig.1.10 

Mi.i 



cr 



-^V-H/+H-5 1 ^+H/-5 1 {^o(^> m 5i ) m U')[fl , V' i? A + ^ A]} (CIO) 



-9vw+s 1 s 2 9h-s 1 s 2 

{B (m 2 H± ,m 2 s ,m 2 H± )[gf A, + gf -4 2 ]}(C.ll) 



v 
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Fig.1.11 



Mi. 



a 



11 



9zw+H-s 1 9szz{B (s,m 2 Sl ,m 2 z )[gf At + gf A 2 }} (C.12) 



where the couplings are given in the table 



(Si, S2) 


9zw+SiS 2 




9h-SiS 2 


(A, H~) 


C/3q 

2c w 


Cf3a 

2sw 


9hH~H- 


(H,H-) 


2cw 


2s w 


9hh~h~ 


(h,G+) 


2c w 


Sf3a 

2s w 


9hH~G+ 


(H, G+) 


c /3q 
2cw 


X0n 
2sw 


9HH-G+ 



Box diagrams 

Fig.1.12 



M Box 



a 



4 s 2 



c /3a s / 3 a m l y{[C (mg,mg, s,m^, 0,m^) + (m 2 , - u)D + (m 2 ^ - u)D 1 



+ (m 2 w -s- t)D 3 ]A 2 - 4DiA} 



(C.13) 



where all the Di functions have as arguments: (m 2 ^, ml, ml, m 2 H± , s, m\, m^y, 0, m 2 ^). The 
corresponding box graph with a heavy neutral scalar is obtained by the substitution 



-M* ox (m h ^m H ) 



(C.14) 



Counter-term diagrams 

Denoting by k^ YiHw{k 2 ) the bare H^W^ 2-point function, Fig. 1.17 - Fig. 1.20, the renor- 
malized one is obtained by adding the counter-term ( |4.2| ). The renormalized non-diagonal 
self-energy T> HW (k 2 ) is thus given by: 



According to the condition 
A is determined by: 



E HW (k 2 ) = E HW (k 2 )-A. (C.15) 
requiring that the renormalized W + -H~ mixing vanishes, 



A = ReT, HW (m 2 H ±) . 



(C.16) 



The self-energy T l Hw( m H±) = ^hw has the following explicit form: 



■spferrn . -ybos 



HW 
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The fermionic part is the sum over the individual fermion families, each one contributing 

= ~ N c^ (m d Y u R d B (m^ + ,mlml) + (m u Y u L d + m d Y u R d )B l (mj I+ ,mlml)), 

(C.17) 

and the bosonic part is given by 

S w = ■^^(ci3 a 9hH-H-swB (m 2 H± ,m 2 h ,m 2 H± ) - c l3a m w s l3a B (m 2 H± ,m 2 h ,m 2 v ) + 
Qirs w 

9hH-G+Sp a SwBQ(m 2 H± ,m 2 h ,m 2 w ) - g H H+H-Sp a s w B Q (m 2 H± ,m 2 H ,m 2 H± ) + 
cp a m w sp a BQ(m 2 H± ,m 2 H ,m 2 w ) + c l3a g HH - G +s w B (m 2 H± ,m 2 H ,m 2 v ) + 
2c fSa g hH + H -s w B 1 (m 2 H± ,m 2 h ,m 2 H± ) + cp a m w sp a B x {m 2 H± ,m\,m 2 w ) + 
ZghH+G+spaSwB^m^mlimlv) - 2g HH + H -sp a s w B 1 (m 2 H ±,m 2 H ,m 2 H ±) - 
c^m w s &a B x {vr? H± ,m 2 H ,m 2 w ) + 2c l3a g HH+G +s w B 1 (m 2 H± ,m 2 H ,m 2 v )) . (C.18) 

Then the amplitudes containing the counter terms, Fig. 1.28, read: 



Ms { vwh) = ^^Re(E HW )[gf A 1 + gf A 2 ] (C.19) 
with gz = —sw/cw, g-y = _ 1 an d gy" ,R defined in (C.l) 

Amplitudes with non-diagonal self-energies 

t-channel W + —H~ mixing 

The generic diagram for this topology is depicted in Fig. 1.13. The amplitude contains only 
the invariant A 2 '■ 



e 2 Im(E^) + Im(E>' 



'OS 

M L13 = __ -v-^ , ■ -y-HwJ A2 (c _ 20) 



m 2 H± - mly 



s-channel W + —H mixing 



The generic diagram for this topology is depicted in Fig.1.15. The amplitude can be 
projected on Ai <2 as follows: 



M L15 = p^-Am 2 w - s) M^Y)+M S w) ( ^ + eLA2) (C 21) 



where gzww = cw/ s w and g^ww = 1 
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s-channel G + —H mixing 

The non-diagonal Goldstone-if ± self-energy is Yiqh evaluated in the same way as done 
before for the W-H ± case {k 2 = m 2 H ±): 

GH - ^GH + ^GH > 

where the fermionic part is obtained by summing over the families with 

^gh = , , {(m u Y u L d - m d Y u R d )A (m 2 d ) + m u (-m 2 d + m 2 u )Y u L d B (m 2 H± , m 2 d , m 2 u ) 

+m 2 H± (m u Y u L d - m d Y u R d )B 1 (m 2 H± ,mlm 2 d )} , (C.22) 
and the bosonic part is given by 

^gh = ^{9hH+H-9hG+H-B (rn 2 H± ,m 2 h ,m 2 H± ) + g H H+H-gHG+H-B {m 2 H± ,m 2 H ,m 2 H± ) 
+9hG+H-9hG+G-B (m 2 H± ,m 2 h ,m 2 v ) + g H G+H-gHG+G-Bo{m 2 H ±,m 2 H ,m 2 w ) - 
2 {A (m w ) + {m h + m H± )B (m H± ,m h: m w ) - 2m H ±B 1 (m H ±, m h , m w ) + 

-(Ao(m^) + (m 2 H + m 2 H ±)B (m 2 H ±,m 2 H ,m 2 v ) - 2m 2 H ±Bi(m 2 H ±, m 2 H , m^))} . (C.23) 

Then the amplitude of the diagram Fig. 1.14 can be written in the form 

ML14 " ( S -<) m 2 H± - m 2 w {dv Al + 9v Ai) ^ 

where gzw+G- — —iriwsw/cw and g^w+G- — n^w- 
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Figure Captions: 



Fig. 1: Feynman diagrams relevant for the e + e~ — > W + H~. Fig. 1.1 — > Fig. 1.11 vertex 
diagrams, Fig. 1.12 box diagram, Fig. 1.13 — > Fig. 1.16 typical diagrams contributing 
the self energies mixing of W + -H~ and G + -H~, self-energies of the mixing W + -H~ 
Fig.1.17 -> Fig.1.20, self-energies of the mixing G + -H~ Fig.1.21 -> Fig.1.25. Fig.1.26 
and Fig.1.27 are the counter term for the mixing G + -H~ and W + -H~ . Fig. 1.28 
are the counter-terms for the photon-IF + -i/~ and Z-W + -H~ vertices. 

Fig. 2: Top-bottom contribution to the integrated cross section as a function of tan/3 in 
Model type I and II for m H ± = 220GeV, and = 500GeV. 

Fig. 3: Top-bottom contribution to the integrated cross section as a function of y/s for four 
values ofm H ± = UOGeV, 185GeV, 300GeV, and 400GeV. Fig. 3.a tan/3 = 0.6 
and Fig. 3.b tan/3 = 1.6 

Fig. 4: Integrated total cross section as a function of y/s for four values of A 3 . 

Fig. 4. a Ci case: m H ± = 220, m H = 180, m h = 90, m A = 80, tana = 1.4, 
tan (3 = 3.6 

Fig. 4.b C 2 case: m H ± = 300, m H = 280, m h = 120, m A = 220, tana = 2.4, 
tan/3 = 1.6 

Fig. 4.c C 3 case: m H ± = 400, m H = 380, m h = 120, m A = 370, tana = 3, tan/3 = 2 
(where all masses are in GeV) 

Fig. 4.d Box contributions to the total cross section as function of the center of mass 
energy for 61,2,3 configurations. 

Fig. 5: Integrated total cross section as a function of tan/3 in the case of the 61,2,3 configu- 
rations with A 3 = 0.1 and for ^/Ts = 500 GeV (Fig5.a), 1 TeV (Fig.5.b) and 1.5 TeV 
(Fig.5.c) respectively. 

Fig. 6: Integrated total cross section as a function of for ^Js = 500, 1000, and 1500 
GeV, with tan/3 = 2, tana = 3, A 3 = 0.1, m H = 180 GeV, m h = 90GeV and 
m A = 80 GeV 
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Fig 4. c (C 3 case ) 
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